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The Wigner-function formalism is a widely used approach to model charge transport in semi- 
conductor nanodevices. Primary goal of the present Paper is to point out and explain intrinsic 
limitations of the conventional quantum-device modeling based on such Wigner-function paradigm, 
providing a definite answer to open questions related to the appHcation of the conventional spa- 
tial boundary condition scheme to the Wigner transport equation. Our analysis shows that (i) in 
the presence of a symmetric device potential the solution of the Wigner equation (compatible with 
given boundary conditions) is not unique, and (ii) also when the solution is unique, the latter is not 
necessarily a physical Wigner function, i.e., the Weyl- Wigner transform of a single-particle density 
matrix. 



I. INTRODUCTION 

Current technology of micro/nanoelectronics pushes 
device dimensions toward space- and time-scales where 
the traditional semiclassical or Boltzmann pictur^32] can 
no longer be a ppli ed, and genuine quantum approaches 
are imperative.'^EI However, in spite of the quantum- 
mechanical nature of electron and photon dynamics in 
the core region of typical solid-state nanodevices -e.g., 
superlattice^ and quantum-dot structures^- the overall 
behavior of such quantum systems is often governed by a 
highly non-trivial interplay between phase coherence and 
dissipation/dephasingjEi the latter being also strongly 
influenced by the presence of space boundaries.^ 
A widely used theoretical tool to account for such 
interplay in semiconductors is the single-particl e de nsity 
matrix operator p for the electron subsystem.!^ The 
time evolution of p is given by the density-matrix 
equation, which involves both a coherent-dynamics 
term and a scattering superoperator encoding the 
energy-dissipative/decoherent interaction mechanisms 
that electrons experience within the host material. 
The density-matr ix a pproach applies to a variety of 
physical problems ,^223 ranging from quantum-transport 
phenomena to ultrafast electro-optical processes in 
"extended systems", i.e., systems extending over the 
whole coordinate space. 

However, such approach cannot be straightforwardly 
applied to nanostructured devices. Indeed, a typical 
nanodevic^ is a "localized system", i.e., a portion of 
material characterized by a well defined volume and by 
spatial boundaries acting as electric contacts to external 
charge reservoirs, as sketched in Fig. [T] Here, z denotes 
the transport direction, I is the longitudinal length of the 
device, the electric contacts being located at z = —1/2 
and z = +1/2. The modeling of a nanostructure device 
thus represents an intrinsically space-dependent prob- 
lem, so that a real-space quantum treatment accounting 
for the presence of quite different spatial regions be- 
comes mandatory. To this purpose, the Wigner-function 
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FIG. 1. (Color online) Schematic representation of a typical 
semiconductor-based quantum device as an open system con- 
nected to two external charge reservoirs. Here, the distance 
between the interfaces is I, and z is the longitudinal transport 
direction. 



formalisrrPE^l ig a, widely adopted framework. Within 
this formalism, the statistical quantum state of the 
electronic subsystem is fully described in terms of the 
Wigner function /(r,k), defined over the conventional 
phase-space (r, k) as the Weyl- Wigner transform of the 
single-particle density-matrix operator p^^ 



Based on the Wigner-function formalism, various ap- 
proaches for the study of quantum-transport phenomena 
in semiconductor nanomaterials and nanodevices have 
been proposed.E^HSl] Qn the one hand, starting from the 
pioneering work by Frensley,!^ quantum-transport simu- 
lations based on a direct numerical solution of the Wigner 
equation have been performed via finite-difference ap- 
proacheJi^, by imposing on the Wigner function the stan- 
dard U boundary condition scheme (see Fig.|2]in Sec. Ill I. 
On the other hand, a generalization to systems with open 
boundaries of t he sem iconductor Bloch equations has also 
been proposed! '^^ * '^^ ! In addition to these two alternative 
simulation strategies -both based on effective treatments 
of relevant interaction mechanisms- Jacoboni and co- 
workers have proposed a fully quantum-mechanical sim- 
ulation scheme for the study of electron-phonon inter- 
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action based on the "Wigner paths" f^H this approach is 
intrinsicaUy able to overcome the standard approxima- 
tions of conventional quantum-transport models, namely 
the Markov approximation and the completed-collision 
limit however, due to the huge amount of computation 
required, its applicability is often limited to short time- 
scales and extremely simplified situations. 

Motivated by a few unphysical result^ obtained 
via the generalized semiconductor Bloch equations 
mentioned above, a recent studj^^ has shown that 
the application of the conventional space boundary 
condition scheme to the Wigner transport equation may 
lead to partially negative charge probability densities, 
unambiguous proof of the failure of such classical-like 
Wigner-function treatment. 

Primary goal of the present Paper is to point out that 
the conventional boundary condition scheme adopted for 
the solution of the Wigner equation exhibits some intrin- 
sic limitations, whose impact may lead to totally unphys- 
ical results, especially in the coherent-transport regime. 
In order to illustrate this aspect, we shall mainly discuss 
a simple but physically relevant one-dimensional system, 
namely a delta-like potential profile. Our detailed analy- 
sis will show that (i) in the presence of a symmetric device 
potential the solution of the Wigner equation (compati- 
ble with given boundary conditions) is not unique, and 
(ii) also when the solution is unique, the latter is not nec- 
essarily a Wigner function, i.e., a Weyl- Wigner transform 
of a single-particle density matrix. 

The Paper is organized as follows: In Sec. |ll] we shall 
summarize the fundamentals of quantum-device model- 
ing, with a special focus on the problem of quantum sys- 
tems with open space boundaries, corresponding, e.g., to 
the case of a semiconductor nanodevice inserted into an 
electric circuit. Section [111] is devoted to the investigation 
of the so-called coherent limit and of the corresponding 
Wigner equation, and will address the main topic of the 
Paper, i.e., the physical versus non-physical nature of the 
solutions of the Wigner equation corresponding to given 
spatial boundary conditions. In Sec. |lV|we shall discuss 
the effect of including energy-dissipation and decoherence 
phenomena. Finally, in Sec. |V] we shall summarize and 
draw a few conclusions. 



tr{W^(r,k)p} 



(1) 



corresponding to the quantum-plus-statistical average of 
the Wigner operatoi^^^ 
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Within such Wigner-function representation the average 
values of charge and current densities at location r are 
given by 



n(r) 



dk 
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and 



J(r) 



dk 

(2^ 



v(k)/(r,k) 
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where v(k) = hk/m* is the velocity of a conduction 
electron with effective mass m* . 

The equation for the Wigner function can be de- 
rived from the equation of motion for the density-matrix 
operatorP 



dp _ 1 
dt ih 



+ r(/5) 



(5) 



Here, the first contribution on the r.h.s. describes the co- 
herent dynamics dictated by a non-interacting Hamilto- 
nian H° , including elastic single-electron scattering pro- 
cesses, while the second term is a linear superopera- 
tor r encoding the energy-dissipative/decoherent scat- 
tering mechanisms that electrons experience within the 
host material. 

By applying the Weyl- Wigner transform ([T]), together 
with its inverse 
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P 



(27r)3 



dr / dk#(r,k)/(r,k) 



(6) 



to the density-matrix equation ([5| , one obtains the equa- 
tion of motion for the Wigner function 



a/(r,k) _ a/(r,k) 



dt 



dt 



9/(r,k) 



dt 



(7) 



II. QUANTUM-DEVICE MODELING BASED 
ON THE WIGNER-FUNCTION FORMALISM 

In order to account for the space-dependent charac- 
ter of a quantum device, a widely employed strategy 
is the Wigner-function treatment of the problem.!^ The 
Wigner function /(r,k) associated to a single-particle 
density-matrix operator p is defined as its Weyl- Wigner 
transforrrfl^l 



/(r,k)= J dr 



r — 



with 



and 



5/(r,k) 



dt 



g/(r,k) 
dt 

where 
e(r,k;r',k') 



J dr' dk'e(r,k;r',k')/(r',k') 



(8) 



= / dr'dk'r(r,k;r',k')/(r',k') , (9) 



(27r)3?i 
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and 

r(r, k; r', k') = tr {^^(r, k) T {w{r', k')) } (11) 

are the single-particle and the scattering superoperators 
written in the (r, k) Wigner picture, respectively. 

For any given set of basis states \a), the single-particle 
density-matrix operator can be expressed in terms of en- 

tri6S Pai ct2 



^ \ai)pa,a^{a2\ 
aiQ2 



(12) 



whereas the Wigner function in ([T]) can be written as 



/(i">k) = H^„2ai(r,k)pa,a2 



(13) 



with 



VFasai (r, k) = y dr'<j>a, + ^ 



^ — 2k-] 



r — 



(14) 

where 0Q(r) = (r|a) denotes the real-space wavefunction 
corresponding to the basis state \a). 

In particular, by choosing as basis states |a) the eigen- 
states of H° corresponding to the energy spectrum e^, 
the noninteracting Hamiltonian can be rewritten as 



H° 



^|a)e„(a| , (15) 

a 

and the density-matrix equation ^ is given by 



dt 



ih 



Paic 



^ ^ ^ aia2,a[a'2Pa[a'2 ■ (16) 



Such set of coupled equations of motion for the density- 
matrix elements Paia2 cire usually referred to as the 
semiconductor Bloch equationsP 

We emphasize that in general the Wigner equation ([t]) 
is non-local both in r and in k. As a consequence, the 
conventional boundary condition scheme adopted to solve 
the semiclassical Boltzmann equation cannot be applied."^ 
However, in order to simplify the problem, the rigorous 
expressions of the single-particle and scattering super- 
operators in Eqs. (Sj) and ^ are often replaced by ef- 
fective/phenomenological operators. In particular, as we 
shall discuss in detail in Sec. |III[ within the conventional 
effective-mass and envelope-function approximations, the 
single-particle operator ([s]) turns out to be local. Further- 
more, introducing a generalized relaxation-time approxi- 
mation (see also SecJIVl), the fully quantum- mechanical 
scattering operator M acquires the fully local fornPi 



a/(r,k) 



dt 



/(r,k)-/°(r,k) 
T(r,k) 



(17) 



describing the effect of dissipation/decoherence, toward 
the thermal-equilibrium Wigner function /°(r,k), via 
a relaxation time T(r,k). It is however important to 
stress that the generalized relaxation-time approxi- 
mation (17 1 does not preserve the basic properties 
of the Wigner function. In particular, a partially 
negative charge probability density [see Eq. (|3])] may 
be obtained if r and/or f° are taken to be r- and/or 
k-dependent. In fact, this problem is already present 
at the level of the semiconductor Bloch equations ( [l6| . 
Indeed, since the Wigner equation ^ is obtained 
via the Weyl- Wigner transform of the density-matrix 
equation, the reliability of the Wigner-equation ap- 
proach depends on the degree of accuracy of such 
density-matrix formalism, which is in turn intimately 
related to the choice of the scattering superoperator F 
in ([5|. Oversimplified approaches accounting for F in 
a phenomenological way lead to the violation of the 
positive definite character of the density-matrix and 
therefore to unphysical conclusions. In contrast, a recent 
microscopic derivatiorP^ of F in terms of Lindblad-like 
scattering superoperator^^ for arbitrary single-particle 
interaction mechanisms ensures the positive-definite 
character of the density-matrix operator p, and enables 
to overcome this serious limitation. Importantly, since 
the Weyl- Wigner transform ([l]) preserves the positive 
character of p, the equation of motion ([t]) obtained 
from a positive-definite scattering superoperator neces- 
sarily preserves all the basic properties of any Wigner 
function, and in particular the fact that the correspond- 
ing spatial carrier density is positive-definite at any time. 



III. THE COHERENT LIMIT AND THE 
WIGNER TRANSPORT EQUATION 

In order to investigate the role played by spatial 
boundary conditions applied to the Wigner function, 
we shall first focus on a fully coherent system/device, 
where energy-dissipation/decoherence processes occur 
over timescales that are much longer than the typical 
timescales induced by H° . In this regime, the density- 
matrix equation ([s]) reduces to the Liouville-von Neu- 
mann equation 



dp _ 1 
dt ih 



H° 



(18) 



For the purpose of the present Paper, it is sufficient 
to consider a one-dimensional system (r,k z,k) de- 
scribed by the envelope-function HamiltoniarP 



H° = K{k) + V{z) 



(19) 



where z and k denote, respectively, the quantum- 
mechanical operators associated to the electronic coordi- 
nate (z) and to the electronic momentum/wavevector (fc); 
generalizations to a fully three-dimensional problem are 
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straightforward. According to the usual prescription of 
the envelope-function theory, the function K in Eq. ( 19 ) 



describes the bulk electronic band, while V describes the 



nanostructure potential profile. The Hamiltonian ( 19 ) 
leads Eq. ( 18 ) to acquire the form 



dp dp 
dt ^ ~dt 



with 



and 



dp 

dt ^ 



dp 
~di 



K 

1 

ih, 



1 

ih 



dp 
dt 



K{k),p 



V{z),p 



(20) 



(21) 



(22) 



Applying the Weyl-Wigner transform to the density- 
matrix equation (20), one obtains the Wigner- function 
equation for f{z,k). In doing that, the Wigner func- 
tion ([T]) can be expressed in two different and equivalent 
ways, corresponding to the momentum (fc) and coordi- 
nate (z) representations, respectively. By setting a = k 
as well as a = z in Eq. ( 13 ) one obtains 



dk'e"''' p[k 



k' 



,fc 



z 



(23) 

\ ^ ^ / 

rf^'e-''^'p(^+^,^-^) • (24) 

These two expressions turn out to be both useful because 



the kinetic and potential contributions (21 1 and (22 1 are 
diagonal in the momentum (fc) and coordinate (z) repre- 
sentations, respectively, i.e. 



dp(kx,ki) 



and 



dt 



dp{zi,Z2) 



K 



K(ki)-K{k2) 
ih 



p{ki,k2) 



dt 



V{Z^) - V{Z2) 



(25) 



(26) 



i) By applying the Weyl-Wigner transform ( 23 1 , as well 
as its inverse. 



= dz- 



27r 



■f{z,k), (27) 



to the kinetic contribution (251, one obtains 

df{z,k) 



dt 



K 



= - / dz' IC{z^ z',k)f{z\k) (28) 



with 



/C(z",/c) = : 



dk'- 



2tt 



K[k + 



-K\k- 



(29) 



The kinetic operator in the Wigner picture, appearing on 
the r.h.s. of Eq. ( p8| , is always local in k and, in general, 
is non-local in z. In particular, by adopting the usual 
effective-mass approximation, 



2 1,2 



h'k 

2m* 



K{k) 

the non-local kinetic operator ( |28| ) reduces to 

df{z,k) 



df{z,k) 
dt 



= -v{k) 



K 



dz 



(30) 



(31) 



where v{k) — hk/m* denotes the effective-mass carrier 
group velocity. Notably, within the effective-mass 



approximation ( 30 ) the kinetic contribution coincides 
with its semiclassical counterpart, i.e., it reduces to the 
standard diffusion term of the Boltzmann equation. 



ii) By applying the Weyl-Wigner transform (24), as 
well as its inverse. 



z z 

z ~\ , z 

2 ' 2 



Az' k 



dk- 



27r 



(32) 



to the potential contribution in ( 25 ) , one obtains 
df{z,k) 



dt 



= - dk'V{z,k-k')f{z,k') (33) 



with 



f -ik"z' 








dz' 


v(^z^ 






J 2tt 









v(z,n = 

(34) 

Oppositely to the kinetic one, the potential operator ap- 
pearing on the r.h.s. of Eq.(33) is always local in z and. 



in general, is non-local in k. For the particular case of a 
quadratic potential 



F(z) = -az-" 



bz + c 



corresponding to the classical force 
dV{z 



F{z) 



dz 



= —{az + 6), 



(35) 



(36) 



the non-local potential operator ( 33 ) simply reduces to 
df{z,k) 



dt 



F{z) df{z,k) 
h dk 



(37) 



Thus, for the particular case of the quadratic potential 



profile ( 35 ) , the potential contribution coincides with its 
semiclassical counterpart, i.e., it reduces to the standard 
drift term of the Boltzmann equation; it follows that 
the non-local character of the generic potential super- 
operator in ( 33 ) vanishes in the presence of a parabolic 
potential only. 
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Equation (38), also referred to as the Wigner trans- 



FIG. 2. The conventional inflow or U boundary condi- 
tion scheme adopted in semiclassical device modeling-, for 
a one-dimensional problem. The value of the Wigner func- 
tion f{z,k) is specified at the boundaries z''{k) of the active 
region, i.e., f{-l/2,k > 0) and f{+l/2,k < 0) are fixed by 
the incoming/infiowing carrier distribution function. 



At this point a few comments are in order. The analy- 
sis performed so far has shown a strongly symmetric role 
between real-space (z) and momentum {k) coordinates; 
this is confirmed by the fact that the corresponding 
equations of motion (each one written within the related 
representation) display the very same mathematical 
structure [see Eqs. (|25]) and (26l]. Moreover, for a phys- 
ical system characterized by an effective Hamiltonian 
quadratic in both the coordinate and the momentum, 
the equation of motion of the Wigner function coincides 
with its semiclassical (Boltzmann) counterpart, thus 
showing the intimate link between the Wigner function 
and the semiclassical distribution. This can also be 
regarded as a formal proof of the fact that, for a particle 
subjected to a quadratic potential, its classical and 
quantum equations of motion coincide, a fundamental 
result originally pointed out by Richard P. Feynman via 
his "path integral" formulation of quantum mechanics.'*^ 



For the microscopic modeling of semiconductor quan- 
tum devices, the effective-mass approximation (30) is 



widely employed, and constitutes a good starting point 
for the description of the bulk band structure. In con- 
trast, for a generic optoelectronic device, the effective 
potential profile V{z) is usually far from the quadratic 
form in (35). As a consequence, within this approxima- 
tion scheme the single-particle superoperator e in ([8| is 
always local in z, and the total (i.e., kinetic plus poten- 
tial) equation of motion for /(z, k) -obtained combining 
Eqs. (31) and (33)- contains a non-local term in k in- 



duced by the potential profile 



a/(z,fc) ,,.df{z,k) 

+ v[k) 



dt 



dz 



+ / dk'V{z,k^ 



-k')f{z,k')^Q . 

(38) 



port equation, describes the time evolution of the one- 
dimensional Wigner function in the absence of energy- 
dissipation/dccohcrence processes. In steady-state con- 
ditions (0/(z, k)/dt = 0) it reduces to 



v{k) 



df{z,k) 
dz 



dk'V{z,k~k')f{z,k') 



(39) 



From a mathematical point of view, in terms of the vari- 
able z, the Wigner Equation (39) is a first-order differ- 



ential equation. In this respect, it is thus similar to the 
(steady-state) semiclassical Boltzmann equatiori^ 



v{k) 



df{z,k) _ F{z)df{z,k) 



dz 



dk 



(40) 



Based on this analogy, outlined in the pioneering work 
by Frensley,^ several quantum-transport problems have 
been treated by following a semiclassical approach, i.e. 
by applying to the Wigner transport equation (39) the 



strat egy commonly adopted for the Boltzmann Equa- 
tion ( 40 ) . Indeed most of these studieJ^ are based on 



a numerical solution of Eq. ( 39 ) , often supplemented by 



an additional relaxation-time term, where one imposes 
on /(z, k) the U spatial boundary condition scheme of 
the semiclassical device modeling. The latter, depicted 
in Fig. [2j consists in requiring that the inflowing Wigner 
function acquires some fixed values at the two contacts 
z = ±V2i and that these values are determined by the 
distribution of carriers incoming from the two reservoirs. 
Explicitly, the values f{—l/2,k) are specified for carriers 
incoming from the left reservoir (fc > 0) and the val- 
ues f{+l/2,k) are specified for electrons incoming from 
the right reservoir (fc < 0). In a compact notation, in- 
troducing z^{k) — sign(fc)^/2 (where sign denotes the 
sign function) , the Wigner transport equation ( 39 ) is thus 
equipped with the fc-dependent boundary condition 



f (fc)^/(z''(fc),fc) 



(41) 



In spite of the numerical efficiency and the seemingly 
straightforward physical interpretation of such semiclas- 
sical approach, it has been recently shown that the latter 
may lead to unphysical results. More specifically, as orig- 
inally pointed out in Rcf. |37| , the application of the con- 
ventional U boundary condition scheme to the Wigner- 
function equation may lead to negative carrier proba- 
bility densities, a clear hallmark that the classical-like 
Wigner-function treatment is not always reliable. In the 
remainder of this Section we shall highlight and discuss 
the conditions for the failure of such a spatial boundary 
condition approach. 

Here we wish to point out that, in spite of the analo- 
gies mentioned above, an important difference emerges 
between the Wigner equation ( 39 ) and the semiclassical 



Boltzmann equation (40). While the latter is local in k, 
the former is not. Indeed, because of the non-local char- 
acter (in k) of the potential superoperator V appearing 
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in (39), the differential equation for one value of k is in 



fact coupled to the differential equations for all other k 
values. The non-locality of V therefore makes the Wigner 
problem intrinsically different from the Boltzmann one. 
We shall show that this mathematical aspect does have 
important physical implications. 



Symmetric potentials: non-uniqueness of the 
solution 



In this subsection we shall analyze in detail the Wigner 
problem, i.e., the steady-state first-order differential 



equation ( 39 ) , equipped with the fc-dependent boundary 



condition (41 1. More specifically, we shall show that, for 



a symmetric potential profile, 



Vi-z) = V{z) , 



(42) 



the Wigner problem may exhibit non-unique solu- 
tions. We shall also specify the conditions for the 
non-uniqueness to occur, whereas in the next subsection 
some explicit examples for such anomalous feature in 
the presence of customary symmetric potentials will be 
provided. 

In order to explain such anomalous non-uniqueness of 
the solution, we start by investigating the general symme- 
try properties of the Wigner transport equation ( |39l). I n 
the presence of a spatially symmetric potential [Eq.(42)], 
the Wigner potential ( 34 ) is easily shown to be antisym- 



metric with respect to the space coordinate, i.e., 

V(-z,fc") = -V(z,fc") . (43) 

This basic symmetry property implies that, for given 
a solution f{z,k) of the Wigner equation (39), also 



f{~z, k) is a solution of the same equation and that, 
because of the linear character of the transport equation, 
any linear combination 

g{z,k)^cifiz,k) + C2f{~z,k) (44) 

is also a possible solution of the Wigner equation (|39]) . 



It is thus useful to rewrite the Wigner problem (39) 
(41) in an equivalent integral forrrPS 



vik)fiz,k)=v{k) f{k) 



dz' 

zi'(k) J -oo 



-t-oo 



(45) 

dk'V{z',k~k')f{z,k') 



and to decompose the generic solution /(z, k) into its 
symmetric and antisymmetric parts 



f{z,k)^fs{z,k)+fA{z,k) 



defined as 



fs/A{z,k) = -{f{z,k)±f{~z,k)) 



(46) 



(47) 



By adding and subtracting Eq. ( 45 ) evaluated in z and 



in — z, one derives the following equations for fs and /a. 



+ 00 



v(k) fsiz, k) = v{k) f\k) ~ I dz' I dk' V(z', k ~ k') (fs{z', k') + fA{z', k')) + 

z''(k) J-oo 
r+ao 

dz' I dk' Viz',k - k') {fsiz',k') + fA{z', k')) 



(48) 



v{k)fA{z,k) = - 



dz' 



+ 00 



dk'Viz',k-k') ifsiz',k')+fAiz',k')) 



Exploiting the antisymmetry (43) of V, the equations for fs and Ja reduce to 



v{k)fs{z,k) = v{k) {f{k)- fA{z\k),k)) 



dz' 



■{k) 



+ 00 



dk'V{z',k-k')fs{z',k') 



v{k)fAiz,k) = dz' 
Jo 



+ C30 



(49) 



dk'V{z' ,k-k') Ja{z' ,k') 



r 



Equation (49) shows that, while the antisymmetric part of the Wigner problem always exist because the trivial 



fA{z,k) fulfills a homogeneous equation that does not function fA{z,k) = is always a solution for the homo- 
involve the symmetric one fs{z,k), the inhomogeneous geneous equation. However, if the equation for admits 
equation for the symmetric part depends on the an- non-trivial solutions, it then admits an infinite set of solu- 
tisymmetric one as well, through its boundary value tions, simply obtained by multiplying any non-trivial so- 
fA{z^{k),k). We notice that purely symmetric solutions lution by an arbitrary constant. In turn, this imphes that 
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the solution to the original Wigner problem (451 is not 



unique. Indeed, let us consider a solution f{z^k) of the 
integral Wigner equation corresponding to the boundary 
condition f^{k) = f{z''{k), k), and its symmetric and an- 
tisymmetric parts fs{z,k) and fA{z,k). Then, for any 
arbitrary constant 7, the function 



flizik) =-/fA{z,k) 



(50) 



is also a solution of the homogeneous equation for 
fA{z,k), i.e. 



v{k)rA{z,k) = - 



dz' 



+ 00 

dk'V{z\k~k')fl{z\k') 

-00 

(51) 

Let us then consider the symmetric solution fg{z, k) ful- 
filling the problem 

v{k) r,iz, k) = v{k) {,f{k) - fliz\k),k)) - (52) 

/+00 
dk'V{z',k- k') f^{z',k') 

Then, the new function 

r{z, k) - fl{z, k) + fl{z, k) + f{z, k) (53) 

is also a solution of the original Wigner problem ( |45|) , 
as can be easily seen by summing up Eqs. (51) and (52) 
and by evaluating f{z = z^{k),k) in Eq. (51). In the 



next subsection we shall provide explicit examples of 
various solutions of the same problem compatible with 
the same boundary condition. 

The lack of uniqueness is thus based on the existence of 
non-trivial antisymmetric solutions of the homogeneous 
equation in (49). This, in turn, is related to the non- 



invertibility of the integral operator entering the Wigner 
problem ( 45 ) . More specifically, by adopting a compact 



notation, the Wigner problem (45) can be written as 



Af = /" (54) 



with 

Af{z,k) = f{z,k) + 



(55) 

dz' dk'y^^-^Jiz'^k') 
'(k) J-00 v[k) 

If A is invertible the solution is univocally written as the 
well-known Neumann series 



00 







(56) 



and, as shown in Ref. |37| . for a symmetric potential 
the result of the above Neumann series is always sym- 
metric. However, if A is not invertible, the subspace of 
antisymmetric solutions previously mentioned leads to 
the non-uniqueness of the problem. 



It is worth emphasizing that the lack of invertibility 
strongly relies on the non-local (in k) character of the 
operator A. Indeed, it is easy to show that for a local 
equation 



v{k) fAiz, k) 



dz'V{z,k)fA{z,k) (57) 



the only antisymmetric solution is the trivial one. It 
may thus seem at first that the lack of invertibility is 
a merely mathematical issue. However, such operator 
is in fact non-invertible for quite customary potential 
profiles that are widely used in device modeling, as we 
shall discuss in detail in the next subsection. 

Finally, we would like to point out some relevant con- 
sequences of the non-uniqueness of the Wigner-problem 
solution in the presence of a symmetric potential pro- 
file. We notice that, because purely symmetric solutions 
of ( 49 ) are always admitted, any numerical solution of 



the Wigner Eq. ( 45 ) based on an iterative scheme where 
the initial guess is a z-independent function [including 
the Neumann expansion Eq. ( |56[ )] will necessarily pro- 
duce only symmetric solutions f{z,k) = f{—z, k) o f the 



Wigner equation. Indeed, any iteration of Eq. (49) does 



not exit the symmetric function subspace, and will there- 
fore lead to a symmetric carrier density profile 



n(z) = p{z,z) 



1 

2^ 



dkf{z,k) =n{—z). 



(58) 



Notably, this result would imply that, for a symmetric 
potential profile, the spatial charge density is always 
symmetric, no matter what the shape of the injected 
carrier distribution /''(fc) determined by the reservoirs 
looks like. Such real-space symmetry, although in 
agreement with the numerical results of the generalized 
semiconductor Bloch equations in Ref. [33] as well as of 
the Wigner-equation treatment in Ref. [37], is in fact 
not to be expected and, even more importantly, it is 
incompatible with the analytical results obtained via ele- 
mentary scattering state calculations (see below). As we 
shall see, this is a clear indication that in this case, due 
to the potential symmetry, the solution of the Wigner 
problem is definitely not unique, and antisymmetric 
solutions do have an extremely important physical role. 
Relevant examples will be given here below. 



1. Example: The case of a delta-like potential 

Let us consider the case of the delta-like potential 

V{z)^AS{z) , (59) 

where A denotes the barrier strength parameter. The 
potential profile is symmetric around z — 0. Within the 
envelope-function and effective-mass approximations, the 
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Schrodinger equation reads 



2m* dz^ 



K5{z) 



(z) = e0(z) 



(60) 



and exhibits a continuous set of doubly-degenerate scat- 
tering eigenstates describing, for a given energy e, injec- 
tion from the left side onto the barrier (left-scattering 
states) 



ikz I ^ ikz 



for z > 



(61) 



and injection from the right side onto the barrier (right- 
scattering states) 



1 f g ikz ^ ^_ ^ikz j-Qj. 2 > /„2\ 

for z < 



V L tfc e 



respectively. In Eqs. (61|-(62), k — -\/2m*e/?i > 0, 



lAj:, 



1 



1 + iAfc ' 



(63) 



denote the reflection and transmission amplitudes, re- 
spectively, 



A. 



(64) 



is an effective barrier strength parameter, and the pref- 
actor ensures the normalization of the above scat- 

tering states over the whole system (device-|-reservoirs) 
length Vl.'^^ 

Let us now consider the Wigner function correspond- 
ing to the generic left-scattering state in (61 1. Generally 



speaking, recalling that the density matrix of a pure state 
in the z-representation is given by 



p{z,z') = (l){z)(l)*{z') , 



the corresponding Wigner function ( 24 ) reduces to 

-ikz' 



f{z,k) 



dz' I 



z 

z+-|e 



z 

z 

2 



(65) 



(66) 



By inserting the explicit form of the scattering state ( |61[ ) 
into the above pure-state prescription, a lengthly but 
straightforward calculation, outlined in the Appendix, 
leads to obtain 



/+(^,fc) 



27r 



2fccos {2{k - k)2 



2nk{k - k) 



Tfe 8{k -k)+ itj,rl \^m{2kz)S{k) 

ei~z)R-k ( sin (2(fc - k)z) sin (2(fc + k)z) 
TT I k — k k + k 



2 cos(2fcz) sin(2fcz) 
k 



(67) 



where Rj. = |rj;p and Tj^ 
transmission coefficients 



are the reflection and 



Ri. 



the barrier-free case. For the case we are presently 
considering -that is, one single scattering state as in 
Eq. (61)- no is simply given by 1/il. 



1 + ^1 



T-, 



1 + ^1 



(68) 



67 1 is non- 



A lengthly but straightforward calculation, summa- 
rized in the Appendix, allows to verify that the Wigner 
function (67) is a solution of the Wigner transport equa- 



It is worth noticing that the Wigner function ( 
symmetric, in accordance with the fact that t 

probabihty density n{z) = l^^p corresponding to the the ddta-like barrier (l59| is now given by 



le electron tion (39 1, where the Wigner potential V(z, k) induced by 



pure left-scattering state (61 1 is given by 



;(z) = no 



1 + 



2Re 



^ — 2ikz 



for z < 
for z > 



4A 

V(z,fc) -TT^ sin(2A:z) 



(70) 



(69) 

and is non symmetric with respect to z. This asymmetry 
has a physically intuitive explanation: since a left 
scattering state describes carrier injection from left, the 
presence of the barrier causes a charge accumulation 
on the left of the barrier with respect to the density 
of the carriers transmitted on the right. Here, as 
well as throughout the whole Paper, Uo denotes the 
(space-independent) charge density corresponding to 



as can be easily verified via its definition in Eq. ( 34 1 . 
At this point, a few comments are in order. On the 



one hand, although the Wigner function (67) does solve 



the Wigner transport equation (39), it is definitely not 
symmetric in z. On the other hand, by applying the 



formal-integration scheme ( 45 1 to the case of the delta- 



like potential ( 59 1 , the resulting Neumann expansion in 



(56) gives always a spatially symmetric solution, neces- 



sarily different from the non-symmetric one in ( 67 ) . It is 



thus clear that, according to the general analysis of the 
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Wigner problem in Sec. Ill A[ for the delta-like potential 
profile considered so far, the solution of the Wigner 



transport equation (39 1 corresponding to given spatial 
boundary conditions is not unique. As an example of 
this non-uniqueness of the solution, Fig. [3] shows three 
carrier-density profiles corresponding to three differ- 
ent solutions of the same Wigner problem, i.e. of the 
Wigner equation with the very same boundary condition. 



2. The subspace of compatible solutions 

As a concrete example of the non-uniqueness charac- 
ter previously discussed, we shall now identify an infinite 
set of solutions of the Wigner equation for the delta-like 
potential compatible with given boundary conditions. In 
order to illustrate this aspect, it is worth replacing the 
(physical) inflow boundary conditions previously intro- 
duced with a simpler (mathematical) scheme where we 
assign/fix the value of the Wigner function in z = (for 
allfc): /''(fc) = /(0,fc). 

Let us now look for particular choices of the linear 



combination ( 44 ) corresponding to vanishing boundary 



values. More specifically, let us consider the solution 



giz,k) = f{z,k) - f{-z,k) , 



(71) 



whose boundary value is always equal to zero, i.e., 
(7(0, k) = 0. It follows that, given a particular solution 
f{z,k) corresponding to the boundary values f^{k) = 
f{0,k), any function 



h{z,k) = f{z,k) + cg{z,k) 



(72) 



is also a solution of the Wigner transport equation (39) 
corresponding to the very same boundary values. This 
confirms that, for any given boundary function f^{k), 
the solution of the first-order differential equation is def- 
initely not unique. 

This result explicitly shows the coexistence of the spa- 
tially non-symmetric solution ( |67[ ) (corresponding to the 
left-scattering state ( [gT] )) and the spatially symmetric so- 
lution provided by the Neumann expansion ( 56 1 (see also 
Fig. [3]). Notice that the spatially symmetric solution 
h{z,k) = h{—z,k) is obtained by choosing c = —1/2. 



Indeed, combining Eqs. (71) and (72), in this particular 
case we obtain 

h{z, k) = \ {f{z, k) + f{-z, k)) = fs{z, k) , (73) 

i.e., the solution h{z, k) is simply given by the symmetric 
part of the original solution /(z, k). 

Let us finally discuss the properties of the spatial car- 
rier density corresponding to the general solution in ( 72 1 ; 



by recalling the relation between the Wigner function and 
carrier density in (58) and denoting with nh{z) the spa- 



tial carrier density corresponding to the solution h{z, fc), 
one obtains 
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FIG. 3. (Color online) Non-uniqueness of the solution of the 
Wig ner problem for the case of a delta-like potential bar- 
rier (59 1. Here, I = 40 nm, ? = 100 meV, and A^. = 1.5, 
corresponding to a transmission coefficient Tj. ~ 0.3. Three 
diflerent space carrier density profiles compatible with the 
same boundary condition, and corresponding to the solution 
set in (741: the scattering-state solution (69 1, correspond- 
ing to c = (black solid curve) , the symmetric part ( 75 I , 
corresponding to c= —1/2 (red dashed curve), and the non- 
physical solution corresponding to c = 1/2 (blue dash-dotted 
curve) (see text). 



As expected, for c — —1/2, one has 

nh{z) = i {n{z)+n{-z)) , (75) 

i.e., the spatial carrier density corresponding to c = —1/2 
is simply given by the symmetric part of the original 
carrier density n{z) (see red dashed curve in Fig. |3|. 
This implies that, by applying the standard integration 
scheme in ( 45 ) to the case of the delta-like potential ( 59 ) 



and using as boundary conditions the ones corresponding 
to the scattering state Wigner function (67), the resulting 
Neumann expansion in ( 56 ) does return a spatial carrier 



density corresponding to the symmetric part of the scat- 
tering state carrier density in ( |69[ ). 

In order to validate this conclusion, we have performed 
a numerical solution of the integral version of the Wigner 
equation (45) via a standard finite-difference technique 



(in terms of a 120 x 120 uniform phase-space grid). Fig- 
ure [4] shows a comparison between our numerical result 
(black solid curve) and the symmetric carrier density 
in ([75| (red dashed curve); as expected, the numerical 
solution (corresponding to the Neumann expansion) 
is symmetric, and -apart from small deviations due 
to discretization as well as to phase-space cut-offs- it 



coincides with the analytical carrier density in ( 75 ) 



At this point, an important comment is in order. 



While the spatial carrier density in (69) corresponds to 



nh{z) — n{z) + c {n{z) — n{—z)) 



(74) 



the physical quantum state (61) and is always positive- 
definite (see black solid curve in Fig. |3]) , the correspond- 
ing carrier density n/j(z) in (74) may be non-physical; 
indeed, it is easy to realize that for large enough values 
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position (nm) 

FIG. 4. (Color online) Case of the delta-like potential barrier 
( 59 1 . The device parameters are the same as in Fig.|3] Com- 
parison between the space carrier density obtained via a nu- 
merical solution of the integral Wigner equation ( |45[) ( black 
solid curve) and the symmetric carrier density in ( |75p (red 
dashed curve). A standard phase-space discretization scheme 
(in terms of a 120 x 120 uniform grid) has been employed (see 
text). 



of the parameter c the latter is no more positive-definite 
(see blue dash-dotted curve in Fig. [s]), essential prereq- 
uisite of any physical quantum state. However, the sym- 
metric carrier density ( 75 ) obtained via the conventional 



Neumann expansion is always positive-definite (see red 
dashed curve in Fig.|3|. 

As we shall see, the presence of such non-physical solu- 
tions of the Wigner transport equation is not necessarily 
ascribed to the non-uniqueness discussed so far. Indeed, 
also for a non-symmetric potential profile -for which the 
solution is always unique (see below)- one may easily ob- 
tain non-physical solutions by imposing arbitrary bound- 
ary conditions according to the conventional scheme of 
the semiclassical theory (see black solid curve in Fig. |8]) ; 
this feature, already pointed out in Ref. [37j, appears to 
be the most severe limitation of conventional Wigner- 
bmction treatments. 



3. Other symmetric potential profiles 

While in the above subsection we have focussed on the 
delta-like potential profile, here we would like to com- 
ment about other relevant cases of (more realistic) sym- 
metric potentials. From its definition in Eq. (34 1, one 



can easily show that, for any symmetric potential profile 
y(z), the Wigner potential reads 



V(z,fc) 



sin(2fcz) V{2k) , 



(76) 



where V is the Fourier transform of V{z). The 



dependence of Eq. ( 76 ) is thus exactly the same as the one 



obtained for the case of a delta-like barrier [see Eq. ( 70 )]. 
In contrast, the k dependence of V(z, k) encodes the in- 
formation about the specific profile shape. In particular. 



for the case of a single rectangular barrier with finite 
width a, 



+ z 



(77) 



with 6 denoting the Heaviside step function, one obtains 

V(z,fc) = -2^V^,sin(2fc.)'^ ■ 

For small wavevectors [ka ^ 1) the potential profile in 
Eq. (78 1 reduces to the delta-barrier case in (70), upon 



identifying V^a — K. In contrast, for /ca > 1, an oscil- 
latory behavior emerges in fc, due to the finite barrier 
width a, which makes the convergence of the convolution 
integral in Eq. ( 39 1 faster. The explicit expression of the 



Wigner function for the single barrier problem can be 
found in Ref. gS]. 

Generally speaking, while for a generic potential pro- 
file the exact Wigner function cannot be computed ana- 
lytically, for any piece- wise-constant potential (e.g., step 
as well as multiple-barrier profiles) the corresponding 
Wigner function can be evaluated analytically by gener- 
alizing the prescription given in the Appendix. However, 
based on the previous remarks about delta-versus-finite 
potential barriers, one can claim that the case of a delta- 
like profile considered in this Paper already contains the 
essential ingredients of the Wigner problem for a sym- 
metric potential. 



B. Non-symmetric potentials: uniqueness of the 
solution 

Let us now discuss the case of a non-symmetric poten- 
tial profile. Based on the symmetry analysis performed 
so far, in this case the Wigner potential is no more an- 
tisymmetric with respect to z, i.e., V(— z, k) ^ — V(z, k). 
This implies that, given a solution f[z,k), the function 
g{z,k) in (711 is not a solution anymore and one is not 



allowed to construct the generic solution h{z,k) in (72 1. 

As a result, for the case of a non-symmetric poten- 
tial the solution f{z, k) of the Wigner transport equa- 
tion (39), corresponding to a given boundary condition 
f^{k), is always unique. Opposite to the symmetric- 
potential case, however, this does not constitute a prob- 
lem/paradox, since in the absence of the real-space anti- 
symmetry in ( 43 ) , the solution obtained via the Neumann 



expansion ( 56 ) is no more spatially symmetric (see black 
solid curve in Fig. [s]). 

In order to test the non-symmetric potential case, we 
consider the following potential profile 



V{z) ^Aidiz)+r]S{z-a)) 



(79) 



With respect to the symmetric delta-like potential in 
(59), Eq. (79) includes a second barrier, displaced by a 



from the first one and characterized by a different barrier 
strength. The relative strength rj represents the asym- 
metry parameter. Also in this case, one deals with a 
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FIG. 5. (Color online) Case of the asymmetric potential pro- 



file (79 1, with a distance a — 5nm between the two delta- like 
barriers, sketched with peaked dashed areas. The device pa- 
rameters are the same as in Fig. [3] The space carrier density 
profile corresponding to the left-scattering state is plotted for 
two different values of the asymmetry parameter -q (i.e. the 
relative barrier height), namely rj = 0.75 (red dashed curve) 
and 77 — 0.25 (blue dash-dotted curve), and compared with 
the ri = case (black solid curve) corresponding to the sym- 
metric delta-like potential in ( |59[ ) (see black solid curve in 
Fig.[3f. 



continuous (and doubly-degenerate) set of propagating 
or scattering states. 

Figure [5] shows the spatial carrier-density profile for 
different values of the relative barrier height rj entering 
the double-delta potential in (79). With respect to the 



case of one single delta-like profile [black solid curve, 
corresponding to the density profile of the scattering 
state (61)], the presence of the second barrier induces a 



coherent-tunneling scenario. Indeed, similarly to the case 
of a resonant-tunneling device (where the transmission 
coefficient displays a resonance peak as a function of the 
incoming electron energy), here -for a given/fixed value 
of e- the tunneling process is dictated/ controlled by 
the relative barrier height 77 [see Eq. ([79|]. For rj = 0.25 
(blue dash-dotted curve) a constructive-interference 
regime emerges, leading to a significant increase of the 
transmission coefficient; in contrast, for 77 = 0.75 (red 
dashed curve) one deals with a destructive-interference 
pattern, leading to a decrease of the tunneling dynamics. 
Notably, for 77 = the potential profile in (79) reduces 
to the symmetric one in ( [59| , and the solution of the 
corresponding Wigner equation is not unique. Thus 
the 77 = carrier density (black solid curve in Fig. [S]) 
-corresponding to the scattering state (61)- is only one 



of the infinite set of solutions in (74) (see also Fig. |3|). 
In contrast, for rj ^ 0,1 the profile is definitely non 
symmetric, and the solution of the corresponding Wigner 
equation (compatible with given boundary conditions) 
is always unique. The related carrier density coincides 
with the squared modulus of its scattering state. 
The relevant conclusion to be drawn from the present 
analysis is that, for the case of a non-symmetric po- 
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FIG. 6. The inflow bou nda ry profile determined by the 
analytical Wigner function ( 67 1 for a delta- like potential bar- 
rier (591, namely f{z — —1/2, k) for k > (left panel), and 



f{z — +1/2, k) for k <0 (right panel). The Wigner function 
is plotted in units of 2tt/^1, and the device parameters are the 
same as in Fig. |4] (see text). 



tential, the conventional boundary condition scheme 
is mathematically consistent. In other words, given 
a Wigner function obtained via the Weyl- Wigner 
transform of a given quantum state [see Eq. (66)], by 



employing the value provided by such Wigner function 
at z^{k) as the boundary conditions f^{k), the conven- 
tional integration scheme returns as unique solution the 
Wigner function itself. 

In spite of such mathematical consistency, the scheme 
may fail to provide physically correct results, as we shall 
see in the next subsection by analyzing the coherent- 
transport regime (see Fig.l8|. 



C. Quantum versus semiclassical treatment of 
spatial boundary conditions 

We would like to discuss now in more detail the role 



played by the boundary conditions (41) imposed on the 
Wigner equation (39). In the first instance we notice 



that, for the case of the delta-barrier profile discussed 
above, the Wigner function may be obtained directly 
from the knowledge of the scattering eigenstates of the 
problem, and it is thus possible to know its boundary val- 
ues in advance. Indeed, as a consistency test, we have nu- 
merically checked [see Fig. [4] that, imposing such bound- 
ary values as boundary conditions for the Wigner equa- 
tion, the obtained solution is indeed the original Wigner 
function. 

At this point, it is vital to stress that the conven- 
tional quantum-device modeling based on the Wigner 
formalism is grounded on a totally different paradigm. 
Since the Wigner function is initially unknown, the spa- 
tial boundary condition f^{k) is treated as an arbitrary 
input function, usually expressed in terms of equilibrium 
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or quasi equilibrium inflowing semiclassical carrier distri- 
butions (see Fig. [2]). As recently pointed out in Ref. [32] . 
such a semiclassical treatment /description of the bound- 
ary function f^{k) seems to be not necessarily compat- 
ible with the quantum-mechanical nature of a genuine 
Wigner function, as confirmed by the highly non-classical 
(i.e., non positive-definite) shape of the boundary condi- 
tions corresponding, e.g., to the scattering state solution 
(67 1. This is clearly shown in Fig. |6j where we report 
the left {k > 0) and right {k < 0) quantum-mechanical 
inflow boundary profile corresponding to the analytical 
Wigner function in (67); as anticipated, opposite to the 
usual semiclassical treatment, here the boundary func- 
tion -corresponding to a left-scattering state of incoming 
wavevector k- involves all k values and, more important, 
is not positive-definite. 

It is however crucial to stress that, thanks to the 
presence of external carrier reservoirs in thermal or 
quasi-thermal equilibrium, the Wigner function of the 
quantum-device electron is expected to be far from a pure 
state. For this reason, in order to better compare the 
rigorous shape of the inflowing Wigner function with the 
usual Fermi-Dirac distribution of the semiclassical the- 
ory (employed in the conventional Wigner-function mod- 
eling) , we consider here the Wigner function correspond- 
ing to a mixed state. To this aim, a thermal average is 
performed as incoherent superposition of the density ma- 



trices |-|-fc)(-t-fc| and \ — k){—k 



and right-scattering states (61 1 and (62), respectively. In 



corresponding to the left- 



more explicit terms, this amounts to assume a single par- 
ticle density-matrix operator of the form 



P 



n 

2^ 



"2^ Jo 



dk\+k){+k\r{-ej,-piL) + 
^dk\-k){-k\r{-e-,-iiR) , (80) 



where /xl and hr denote the chemical potentials of the 
left and right reservoirs, and f° is the corresponding 
quasiequilibrium Fermi functions. By applying to the 
mixed-state density matrix ( [80| the Weyl- Wigner trans- 
form in (24), the corresponding Wigner function comes 
out to be 

fiz.k)^— I dkf+iz,k)nE-,^fiL) + 



'0 

2^ Jo 



dkj^{z,k)nEt,~iiR) 



u) 



where, as discussed in the Appendix, the right-scattering 
Wigner function is obtained from the left-scattering 
one in \Q7\ via the transformation z, fc — >■ — z, —k: 
f^(z,k) = f^{-z,^k). 

In order to quantify the impact of the above thermal 
average (with respect to the pure-state result in Fig. |6]) , 
we have evaluated the inflowing part (fc > 0) of the 
Wigner function (81) on the left boundary (z = —1/2) 
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FIG. 7. Case of the delta-like potential barrier (59 1. The 
value of the thermally averaged Wigner function in (81 1 at 
the left boundary is plotted as a function of the wave vec- 
tor k (dashed curves) and is compared to the semiclassical 
assumption of a Fermi-Dirac distribution (solid curves), for 
three different temperature values: T — 300 K (a), T = 30K 
(b), and T = 3K (c). For all three cases we have assumed a 
chemical potential hl = iksT (see text). 



Figure [7] shows a comparison between the left-contact 
Wigner function in (81) (dashed curves) and the corre- 
sponding Fermi-Dirac distribution (solid curves) at three 
different temperatures: T = 300K (a), T = 30K (b), 
and r = 3K (c); for all three cases we have assumed a 
chemical potential fJ^L = 'iksT. As one can see, while 
at room temperature [panel (a)] the two curves coincide 
over a large range of k values, for low temperatures 
[panels (b) and (c)] the value of the Wigner function 
on the left boundary signiflcantly differs from the 
Fermi-Dirac distribution, unambiguous proof of the 
failure of a classical-like boundary condition treatment 
in the low-temperature limit. 



assuming carrier injection from left only (/i^ — > —00) 



The above remarks about quantum-versus-classical 
boundary values lead to a crucial question: for a given 
semiclassical-like (i.e., positive-definite) boundary condi- 
tion f^{k), under which condition is the corresponding 
solution f{z,k) of the Wigner transport equation (39) 
physically acceptable? In order to address this impor- 
tant issue, we have performed a numerical solution of 
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FIG. 8. (Color online) Spatial carrier dens ity corresponding 
to the asymmetric potential profile in ( |79[ ) (here, the device 
parameters are the same as in Fig. [5] and the relative bar- 
rier height is = 0.75): comparison between the numerical 
result obtained assuming a monoenergetic injection from the 
left (black solid curve) and the analytical result corresponding 
to its scattering state boundary condition (red dashed curve) 
(here, the numerical result has been obtained via the very 
same phase-space discretization scheme and parameters em- 
ployed in Fig. |4|. Similarly to the result reported in |37| . in 
the monoenergetic-injection case (black solid curve), negative- 
value regions are found, a clear fingerprint of the non-physical 
nature of such Wigner-equation solution (see text). 



the integral Wigner equation (45) for the case of the 



asymmetric potential in (79) replacing the exact (i.e., 
quantum-mechanical) boundary condition with a simple 
(semiclassical) monoenergetic carrier injection from left, 
i.e., 



f{k) cx S{k-k) 



(82) 



The result of such classical-like calculation (black 
solid curve), together with the corresponding quantum- 
mechanical (i.e., scattering state) one (red dashed curve), 
is reported in Fig. |8] As one can see, in addition to the 
presence of an anomalous amount of charge on the right- 
hand side of the device with respect to the scattering 
state result (red dashed curve) ^ the carrier-density pro- 
file corresponding to the monoenergetic injection (solid 
curve) displays a significant negative-value region and 
thus comes out to be definitely non-physical, i.e., no more 
positive-definite. This result is qualitatively similar to 
the one reported in for the case of a cosine-like po- 
tential, thus confirming the strong physical limitations of 
the conventional Wigner-function modeling. 

In order to point out explicitly the intrinsic limita- 
tions of the Wigner-function formalism combined with 
the classical-like boundary scheme (see Fig. [2]), we finally 
focus on the electric current flowing through the quan- 
tum device; the latter may be expressed in terms of the 
Wigner function f{z,k) as 



/(z) oc 



+ 00 



v{k)f{z, k)dk . 



(83) 



Thanks to the charge continuity equation, in steady- 
state conditions the electric current is z-independent 
[I{z) = /o],'^^ and can thus be evaluated in any space 
point. To this end, let us apply Eq. (83) to the left- 
boundary point z — —1/2^ by splitting the integration 
domain into negative and positive k values 



lo oc 



v[k)f 



I 

2'' 



k dk 



+ CX3 



v{k)f 



I 



k]dk. 
(84) 

For the sake of simplicity, let us consider the case of a 
symmetric potential, for which there always exists at 
least one spatially symmetric solution, i.e., /(z, fc) = 
f{—z,k). Recalling that, a,t z — —1/2 for k > and 
at z = +Z/2 for /c < 0, the Wigner function coincides 
with the inflow boundary function f^{k), and using the 
space symmetry of the Wigner function (between —1/2 
and +1/2), Eq. (84) can simply be rewritten as: 



oc 



+ 00 



v{k)f{k)dk . 



(85) 



This equation would imply that, at least for the case of a 
symmetric potential;^ the electric current is determined 
by the boundary values only, and is fully independent of 
the shape of the device potential profile. In particular, 
the value of the current turns out to be the same for a 
potential-free ballistic device as well as for an infinitely 
high potential barrier. This unphysical conclusion ex- 
plicitly shows the intrinsic limitations of the conventional 
Wigner-function treatment based on the classical-like in- 
flow boundary condition scheme. We shall come back to 
this serious drawback in Sec. ITVl 



D. Mathematical versus physical solutions 

From the analysis performed so far, one is forced to 
conclude that the set of mathematical solutions f{z,k) 
corresponding to the set of all possible semiclassical 
(i.e., positive-definite) boundary conditions f^{k) con- 
tains also non-physical solutions, i.e., solutions that do 
not correspond to the Weyl- Wigner transform of a den- 
sity matrix. In order to gain more insight about physi- 
cal versus non-physical solutions, we rewrite the general 
relation ( 13 ) between single-particle density matrix and 



Wigner function for the one-dimensional problem; denot- 
ing with (/'q(z) the single-particle eigcnfunctions of the 
effective Hamiltonian (IGOI, we have 



(86) 



with 



Wa^ai (z, k) ^ I dz'(j)^^ ( ^ + 2" ) ^ '^"2 I ^ 2 



(87) 
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Comparing the explicit form of the above Weyl-Wigner 
unitary transformation with the definition of the pure- 
state Wigner function in (|66[ ), it is easy to reahze that 
the diagonal elements in ( |87[ ) coincide with the Wigner 
function corresponding to the single-particle state a: 

Waa{z,k) = faiz,k). 

In view of the linear character of the Wigner trans- 
port equation, given two Wigner functions fa{z,k) and 
fl3{z, k) corresponding, respectively, to the single-particle 
wavefunctions 4>a{z) and <j)p{z), the linear combination 



/(z, k) = Cafa{z, k) + cpfp{z, k) 



is also a mathematical solution of the same equation, 
whose spatial carrier density [see Eq. ( 58 )] will be of the 
form: 



n{z) = CaUaiz) + Ci3np{z) 



(89) 



In spite of the positive-definite character of the spatial 
charge densities Ua and n^, it is obvious that a proper 
choice of the free coefhcients and cp may give rise to a 
partially negative charge distribution, which corresponds 
to a non-physical solution. On the other hand, the pres- 
ence of given spatial boundary conditions imposes addi- 
tional constraints on the two coefficients Ca and cp, thus 
reducing the set of available solutions. 

Generally speaking, it is a matter of fact that the set of 
solutions of a given differential equation is usually larger 
than the physically acceptable ones. A simple and rele- 
vant example is provided by the steady-state version of 
the Liouville-von Neumann equation (18), i.e.. 



H\p =0 



(90) 



applied to a two-level system. Combining Eqs. (12) and 
(15), it is easy to verify that the two following density 
matrices 



Pa 

















I) 



(91) 



(corresponding to the upper (b) and lower (a) quantum 
state of the system) are both solutions of the above 
Liouville-von Neumann equation. Once again, in view 
of the linearity of the density-matrix equation (90), any 



linear combination of the above pair of solutions is also 
a solution; the simplest example is 



Pa ~ Pb 



1 

-1 



(92) 



However, exactly as for the linear combination (88), the 
last solution is non-physical, since it does not correspond 
to a density matrix; indeed any density matrix is by def- 
inition positive-definite, i.e., its eigenvalues are always 
non- negative (see below). 

In order to better understand the link among the sys- 
tem density matrix Pqiqj, the Wigner function f{z,k), 
and the corresponding boundary function f^{k), let us 



examine in more detail the Weyl-Wigner transform in 
(86). In view of the Hermitian character of the density 



matrix {Pa^^^ ^ Pa2cn)' there it always exists a unitary 
transformation 



U, 



/3a 



(/3k 



such that 



(93) 



(94) 



i.e., able to diagonalize the density matrix; it follows that 
within the new basis 



(95) 



the density matrix is diagonal with positive-definite 
eigenvalues ep > 0. By applying the unitary transfor- 
mation ( 93 ) to the Weyl-Wigner transform in ( 86 ) , one 
has 



f{z,k) = '^epfi3{z,k) 



(96) 



i.e., the corresponding Wigner function can be written as 
a linear combination of the various Wigner functions cor- 
responding to the new quantum-mechanical states /3 in 
(95); in view of the positive- definite nature of the eigen- 



values ep, the above linear combination can be regarded 
as a statistical average (i.e., a so-called mixed state) of 
the various pure-state Wigner functions fp(z,k). This 
result clearly shows that, starting from a density matrix, 
the resulting Wigner function f{z,k) in (86) is always 
physically acceptable. This is easily confirmed by eval- 
uating the spatial carrier density corresponding to the 
Wigner function in (|96[): 



n{z) 



•pnpiz) 



(97) 



Indeed, the latter, being a linear combination of positive- 
definite functions np{z) with positive-definite coefficients 
6/3, is positive-definite as well. 

Let us consider once again the conventional bound- 
ary condition scheme (see Fig. [2]) employed for the sim- 
ulation of quantum devices with open spatial bound- 
aries. In this case the numerical approach proceeds as 
follows: given the device potential profile V{z) we arbi- 
trarily choose/fix the value of an inflowing semiclassical 
(i.e., positive-definite) carrier distribution, regarding the 
latter as the spatial boundary f^{k) = / (z''(fc), fc) of the 
unknown Wigner function /(z, k). 

For any given (real and positive-definite) boundary 
function f^{k), the numerical solution of the integral 
equation ( |45[ ) always provides a real solution f(z,k). 
However, such mathematical solution does not necessar- 
ily have a physical counterpart. In order to provide a 
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definite explanation of tliis crucial point, it is vital to re- 
call that the set of Weyl-Wigner functions WqjQ2(z, /c) 
constitutes a complete set over the phase-space z, k, i.e., 

^ E k)W^^o., (z', k') = 5(z - z')6{k - k') . 

Qia2 

(98) 

This means that, regardless of the particular 
shape/nature of the solution f(z,k), the latter can 
always be expressed as a linear combination of the 
Weyl-Wigner basis functions W^a^^j^ (z, fc): 



(99) 



Since the solution /(z, fc) is always real and since the 
Wigner operator W is Hermitian, the coefficient matrix 
Caia2 is necessarily Hermitian as well: c^_^^^ = c*^^^. 
Similarly to the previously considered case of the density 



matrix Paia2 [see Eqs. (93) and (94)], one can always find 
a unitary transformation 



such that 



(100) 
(101) 



i.e., diagonalizing the coefficient matrix Ca-^a^ in (99); it 
follows that within the new basis 



(102) 



the coefficient matrix c is diagonal with eigenvalues e^. 
By applying the unitary transformation ( 100 ) to the lin- 
ear combination (|99|, one obtains 



(103) 



f{z,k) = ^ e^f^{z,k) , 



i.e., the generic solution /(z, k) can be written again as a 
linear combination of the various Wigner functions cor- 
responding to the new quantum-mechanical states /? in 
( 102 1. However, opposite to the density-matrix case pre- 



viously considered, here the eigenvalues ej are not neces- 
sarily positive-definite; this implies that the mathemati- 
cal solution /(z, k) is not necessarily a Wigner function, 
i.e., it may be non-physical. This is easily confirmed by 
evaluating once again the corresponding spatial carrier 
density: 



n{z) = e^n^(z) 



(104) 



Indeed, although this is again a linear combination of 
positive-definite functions, is not necessarily positive- 
definite due to the possible presence of negative eigen- 
values e^. This is exactly the explanation of the non- 
physical result reported in Fig. |8] 



IV. INCLUDING ENERGY-DISSIPATION AND 
DECOHERENCE PHENOMENA 

In this section we shall discuss how to extend the 
coherent-limit treatment considered so far in order 
to account for energy-dissipation as well as decoher- 
ence phenomena induced by non-elastic scattering pro- 
cesses. To this purpose, the simplest model is the well- 
known relaxation-time approximation;'^ within the single- 
particle density-matrix formalism introduced in Sec. [H] 
[see Eq. ([5|] the latter amounts to adopting a scattering 
superoperator 



r (p) = — — , 

T 



(105) 



where 



P = 



^ J dk|k)r(ek-A*°)(k| (106) 



is the equilibrium density-matrix operator (expressed via 
the corresponding Fermi-Dirac distribution /°) and r de- 
notes the relaxation time. By applying to the scatter- 



ing superoperator ( 105 ) the Weyl-Wigner transform in 



([T]), it is easy to obtain the relaxation-time term in (17), 
where the equilibrium Wigner function coincides with the 
Fermi-Dirac distribution f°, i.e., /°(r,k) = /°(ek — ^°), 
and the Wigner-function relaxation time T(r, k) coin- 
cides with the (space and energy independent) relaxation 
time r. The latter can be regarded as a sort of effective 
or average coherence time, and is mainly determined by 
carrier-phonon as well as carrier-carrier scattering. It is 
worth stressing that, opposite to the general relaxation- 
time prescription in (17), the present density-matrix for- 



mulation via a constant relaxation time r ensures the 
positive-definite character of the single-particle density 
matrix and therefore of the corresponding Wigner func- 
tion. 

For the one-dimensional case considered above, the 
relaxation-time approximation implies the appearance of 
an additional contribution to the steady-state Wigner 
equation (39), leading to the generalized Wigner equa- 
tion 



dz J 

fiz,k)-nk) 



(107) 



with f°{k) = f° {ek — fJ-") denoting the equilibrium Fermi- 
Dirac distribution induced by the host material. Oppo- 



site to its coherent version in ( 39 ) , the above transport 
equation does not exhibit the z — > — z symmetry dis- 
cussed in Sec. |III A[ It follows that, also for the case of a 
symmetric potential profile, its solution is always unique. 

In order to study the interplay between coherence 
and dissipation/decoherence, we have investigated carrier 
transport through the delta-barrier potential ( 59 ) in the 



presence of a room-temperature carrier injection from the 
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FIG. 9. (Color online) Spatial carrier density for the case of a 
delta-like potential. Comparison between the result obtained 
by an analytical approach to the problem (dashed curve) and 
via a numerical solution of the Wigner equation (107 1, for 



different values of the relaxation time r (solid curves), in the 
presence of a room-temperature carrier injection only from 
the left {fiR — >■ — oo) for the same device and simulation pa- 
rameters considered in Fig.|4| and for = '^ksT (see text). 



left reservoir only [ji^ — ?> — oo). While in the coherent 
regime, the Wigner function for this problem is known 
analytically [see Eqs. (81 1 a nd ([67| ], a numerical solu- 
tion of the Wigner equation ( 107[ ) has been obtained for 
different values of the relaxation time r. We have chosen 
the same device and simulation parameters considered in 
Fig. |4| and = '^keT as in Fig. [7] 
Figure [9] shovifs the obtained spatial carrier density. The 
dashed curve represents the coherent-transport charge 
density, and has been obtained directly from the Wigner 
function of the thermal state ( 80 1 , built out of scatter- 



ing states (61|-(62). It is therefore immune from the 



non-physical behaviors of the Wigner equation pointed 
out above. Explicitly, the density is given by the ther- 
mal average of the carrier density in ( 69 1 . In this coher- 



ent regime, one recovers the carrier density predicted by 
Landauer-Biittiker theorjSB while on the right hand side 
of the barrier the density is given by the fraction of the 
carriers injected from the left reservoir that is transmit- 
ted across the barrier, on the left hand side the density 
accumulation is determined by the fraction of carriers 
reflected back to the left reservoir. Notice that, witli 



respect to the density profile in Eq. ( 69 ) for a monochro- 



matic injection at fc, here the thermal state (80) leads to 
effective transmission and reflection coefficients Tg/ / and 
Reffi obtained as thermal average of Tj,'s. Besides, the 



oscillatory contribution in Eq. ( 69 ) averages out far from 
z = 0. 

The solid curves in Figure [9] correspond to three different 
values of the relaxation time: r = 500 ps, t = 5 ps, and 
T = 50 fs. These timescales have to be compared to the 
average transit time, which is given by the ratio between 
the device length and the dissipation-free carrier drift ve- 
locity, and is of the order of 100 fs. Thus, for r = 500 ps 
the impact of energy relaxation and decoherence is ex- 



FIG. 10. Charge current / as a function of the relaxation 
time r (in units of its potential- and dissipation-free value 7o) 
(solid curve) compared to the coherent- limit (r — )■ oo, dashed 
curve) current corresponding to the analytical charge density 
predicted by the Landauer-Biittiker theory (see dashed curve 
in Fig. [gjl (see text). 



pected to be definitely negligible, and the coherent limit 
previously considered should be recovered. However, the 
charge density obtained via a numerical solution of the 
Wigner equation (107) turns out to be significantly dif- 
ferent from the dashed-curve one. Here, in spite of the 



presence of the potential barrier, all injected carriers are 
transmitted, and no refiection takes place. This confirms 
the fact that for a symmetric potential the value of the 
current provided by the solution of the Wigner equation 
is totally independent of the shape of the device potential 
[see Eq. (85)]. In contrast, for smaller values of the re- 
laxation time the impact of dissipation and decoherence 
becomes significant: a decrease of t (see solid curves in 
Fig. [9]) leads to a progressive decrease of the current, as 
shown in Fig. 10 (solid curve), which corresponds to an 



effective reflection of the injected carriers back to the left 
reservoir, induced by the relaxation-time term in (107). 

The coherence-versus-dissipation scenario described so 
far is fully conflrmed by the electronic-current analy- 
sis presented in Fig. |10| here, we report the current 
/ as a function of the relaxation time r, in units of 
its potential- and dissipation- free value /o, (solid curve) 
compared to the dissipation-free current predicted by the 
Landauer-Biittiker theory (dashed curve) corresponding 
to the dashed-curve charge-density profile in Fig. [9] As 
anticipated, in the coherent limit (r — > 00), while within 
the Landauer-Biittiker formalism the presence of the po- 
tential barrier leads to a significant attenuation of the 
current (///o — 0.6, dashed curve), the dissipation-free 
current obtained via the Wigner equation (107) coincides 
with its potential- free value lo, i.e. with the unphysi- 
cal result described by Eq. (85). For decreasing values 
of T -corresponding to an increased impact of energy- 
dissipation/decoherence processes- one observes a pro- 
gressive reduction of the current (see solid curve). 

From the numerical analysis reported so far, one con- 
cludes that the conventional Wigner-function treatment 
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leads, in general, to an overestimation of tunneling-like 
phenomena; such overestimation, particularly severe in 
the coherent-transport limit, may be quantitatively miti- 
gated by the presence of non-elastic scattering processes. 
It is however important to point out that, from a funda- 
mental point of view, the conventional Wigner-function 
treatment of coherent transport is intrinsically incom- 
patible with well established results of the Landauer- 
Biittiker formalism!^ 



on the other hand, we shall verify that such Wigner func- 
tion fulfils the corresponding Wigner equation. 

We start by introducing the general prescription for 
the analytical evaluation of the one-dimensional pure- 



state Wigner function in (66). To this end, we shall limit 



ourselves to quantum-mechanical states whose wavefunc- 
tions have different analytical expressions on the left 
(L) and on the right (i?) of the space-coordinate origin 
{z = 0), i.e., 



SUMMARY AND CONCLUSIONS 



(^) 



(t>L{z) for 
(j)R{z) for 



z < 
z > 



(108) 



In this work we have pointed out and explained some 
intrinsic limitations of the conventional quantum-device 
modeling strategy based on the well-known Wigner- 
function formalism. More specifically, we have provided a 
clear and definite answer to a few open questions related 
to the application of the conventional space-boundary 
condition scheme to the Wigner transport equation. Our 
detailed investigation -consisting in a proper combina- 
tion of analytical and numerical results- has shown that 
(i) in the presence of a symmetric device potential pro- 
file, the solution of the Wigner equation (compatible with 
given boundary conditions) is not unique [see Eq. (72)], 



and (ii) also when the solution is unique, the latter is 
not necessarily a physical Wigner function (see Fig. |8]), 
i.e., a Weyl- Wigner transform of a single-particle density 
matrix. 

From a physical point of view, such intrinsic limitations 
of the standard (i.e., semiclassical) boundary condition 
scheme applied to the Wigner transport equation can be 
summarized as follows: The essentially wrong ingredi- 
ent in the conventional treatment is the artificial space 
separation between device active region (|z| < 1/2) and 
external reservoirs (|z| > 1/2) (see Fig. [2]). Indeed, the 
latter is intrinsically incompatible with the well-known 
non-local character of quantum mechanics. 

In order to overcome the basic limitations of the 
Wigner-function modeling discussed so far, the cru- 
cial step could be to replace the local (i.e., classical- 
like) boundary condition-scheme treatment of the device- 
reservoir interaction with a fully non-local approach; 
to this end, in order to ensure/maintain the positive- 
definite character of the electronic density matrix, a 
possible strategy is to describe the system/device- 
environment/reservoir interaction via a Lindblad-like 
coupling term.**^ This task is beyond the purpose of the 
present work, and will be discussed elsewhere.'^ 



ANALYTICAL EVALUATION OF THE WIGNER 
FUNCTION FOR A DELTA-LIKE POTENTIAL 
PROFILE 

The goal of this Appendix is twofold: on the one hand, 
we shall discuss the analytical derivation of the Wigner 
function corresponding to the delta- like potential in (591; 



This applies to any potential profile of the form 

v{z) = K5{z) + VoOiz) , (109) 

which includes, as particular cases, the delta-like poten- 
tial in (59) [see Eq. (61 1] as well as the standard step- 
potcntial (not considered in this work). 

In order to evaluate the explicit form of the Wigner 
function in ([66]) , the key step is to perform the integration 
over z'; to this end, for any given value z the arguments 
of the two wavefunctions may assume negative (left) as 
well as positive (right) values according to the value of 
z' . In particular one has 



z' < ~2z 


z + 


z 
2" 
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z' > -2z 


z + 


z' 
2" 
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z' >2z -> 


z — 
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z' <2z -> 


z — 


z' 
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> 



(110) 



According to the above set of inequalities, the integration 
domain in (66) (— oo < z' < -|-c») needs to be split into 



three different subdomains. More specifically, for z > 
we have 



f{z,k) 



-|-oo 



dz'e 
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while for z < we have 



fiz,k) 



dz'e 



-ikz' 



— 00 

+ 22 
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(112) 



Taking into account that for both cases (z > and z < 0) 
the last integral is exactly the complex conjugate of the 
first one, i.e., 



-2\z\ 



+ ■00 



(113) 



and that 



f{z')dz' = / f{z')dz' - / f{z')dz' , (114) 

2\z\ Jo Jo 



the two results in (111) and (112) can be combined as: 




-Oi-z) 
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(115) 



We stress that the above prescription can be easily ex- 
tended to any piece- wise-constant potential, like, e.g., 
multi-step as well as multi-barrier profiles. 

For the particular case of the delta-like potential pro- 



file (59), the explicit form of the left (z < 0) and right 
{z > Ojpart ((/)i and (pji) of the electron wavefunction is 
provided by the scattering states in ([6T|)-(62). In partic- 
ular, by inserting into Eq. ( 115[ ) the left scattering state 
(61 ), after a lengthly but straightforward calculation one 
obtains the Wigner function in Eq.(67). It is possible 



to show that the Wigner function corresponding to the 
right scattering state ( 62 ) can simply be obtained from 
the left-scattering one in (67 1 by replacing z with — z as 
well as k with —k: f^{z,E) = f^{—z,—k). To this aim. 



it is crucial to notice that the application to Eq. ( 66 ) of 
the Wigner-space transformation z, fc — >■ — z, —k is equiv- 
alent to replacing 4>{z) with 0*(— z), the very same wave- 
function transformation linking left and right scattering 
states. 

As a final step, we verify that the Wigner function 



( 67 ) is a solution of the corresponding Wigner equation. 



To this aim, it is useful to rewrite the result in (67) by 



expressing the various reflection and transmission ampli- 
tudes/coefficients in terms of the dimensionless parame- 
ter Afc [see Eqs. ^ and ([64])] 



/^(z, k) = ^ ^^^^ I^S{k -~k)-X-, sm{2kz)5{k) + 



2Afcfccos {2{k - k)z) 
27rfc(fc - k) 



6i(-z)A| /sin(2(fc- fc)z) sin (2(A: + fc)z) 2 cos(2fcz) sin(2fcz) ^ 



k-k 



k + k 



(116) 



By inserting the potential superoperator ( 70 ) correspond 



ing to the delta-like barrier profile ( 59 ) into Eq. ( 39 1 , the 



explicit form of the Wigner equation comes out to be 
df{z,k) _ 4A 



v{k) 



dz 



2Trh 



dfc'sin(2(fc-fc')z)/(z,/c') . 

(117) 



In order to verify that the Wigner function ( 67 1 is indeed 



a solution of the above Wigner transport equation, let us 
now evaluate separately its kinetic and potential terms. 
As far as the kinetic contribution is concerned, after a 
tedious but straightforward calculation one obtains 



v{k) 



df{z,k) 
dz 



4Xj.hk 



Q m* 



(sin (2(fc - k)z) - 0{-z)\j, (cos (2(fc + k)z) - cos (2(fc - fc)z))) 



(118) 
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Let us now come to the potential contribution in (117 1. Wigner function (67), again after a tedious but straight 



By inserting the expHcit form of the scattering state forward calculation one obtains 



4A 
2^ 



d/fc'sin(2(fc- k')z)f{z,k') = - 



4A 



(sin (2(fc - k)z) - 9{-z)X-k (cos (2(fc + k)z) 



cos {2{k ~ k)z))) . 

(119) 



By inserting the explicit forms of the kinetic and poten- 
tial terms in (118) and (119) into the Wigner transport 



equation (117), we clearly see that the scattering state 



Wigner function (67) is indeed a solution of the Wigner 
transport equation for A^. = m* K/h^k, the very same 



prescription in ( 64 ) obtained via a direct solution of the 



Schrodinger equation. 
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